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Abstract 



New formulas on the inverse problem for the continuous skew- 
self-adjoint Dirac type system are obtained. For the discrete skew- 
self-adjoint Dirac type system the solution of a general type inverse 
spectral problem is also derived in terms of the Weyl functions. The 
description of the Weyl functions on the interval is given. Borg- 
Marchenko type uniqueness theorems are derived for both discrete 
and continuous non-self-adjoint systems too. 



1 Introduction 



Skew-self-adjoint Dirac type system 
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I p is the p x p identity matrix and v is p x p matrix function, is a classical 
object of analysis. It is also called Zakharov-Shabat or AKNS. In particular, 
system (jl.ljl and its discrete analog 

w k+1 (\)-w k (\) = ~c k w k (\), c k = c; = q-\ k = o,i,... (1.3) 

are auxiliary linear systems for many important integrable non-linear equa- 
tions |TJ [7J H2]. Various results and references on the scattering theory for 
system (jl.lj) one can find in [7j. Weyl functions for this system on the in- 
terval and semi-axis have been introduced in pH]- (Weyl functions are also 
called Weyl-Titchmarsh or M-functions in the literature.) The existence and 
uniqueness of the Weyl function for system (jl.lj) with a bounded on the semi- 
axis potential was proved in , and inverse problems in terms of the Weyl 
functions on the interval and semi-axis have been solved. Interesting recent 
spectral results on the non-self-adjoint (especially skew-self-adjoint) Dirac 
type systems one can find in HH Ej . See also further developments of 
the Weyl-Titchmarsh theory with respect to system (jl.lj) in 0^111 EE]- Weyl 
functions are successfully used in solving initial-boundary value problems for 
integrable nonlinear equations (see [21 EH E3 EE El 1221 123)- Moreover, 
Borg-Marchenko type results on the potentials coinciding on a part of the in- 
terval, where the systems are defined, in terms of the asymptotics of the Weyl 
functions of these systems are of great current interest, and the self-adjoint 
case have been considered recently in jH El El CHI EH Ho]. 

In Section 2 we shall obtain Borg-Marchenko type result and some new 
formulas on the inverse problem for skew-self-adjoint system (jl.lj) . 

The largest Section 3 of the paper is dedicated to the discrete case - 
system ()1.3j) . In the case of the discrete system (jl.3j) on the semi-axis explicit 
procedure to recover potential from the rational Weyl function one can find in 
[T2*] . Here we shall solve inverse problem for the discrete case on the interval 
without requiring the Weyl function to be rational. The set of the Weyl 
functions will be completely described in terms of the Taylor coefficients. 
Borg-Marchenko type theorem will be derived too. 

We denote the complex plane by C, upper (lower) semi-plane by C + 
(C_), and we denote the set of bounded operators acting from Hi into H 2 
by {Hi, H 2 }. 
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2 Continuous case 

In this paper we shall consider systems (jl.ljl with locally bounded poten- 
tials on the intervals [0, /] and [0, oo). In other words we assume that the 
inequalities 

\\v{x)\\<M (0<x<0 (2.1) 

are true. Normalize the m x m (m = 2p) fundamental solution u of system 
by the condition u(0, A) = I m . In view of (jl.lj) it is immediate that 

u(x, X)*u(x, A) = u(x, X)u(x, A)* = I m for A = A. (2.2) 

We shall use also the notations for the rows of u(x, 0): 

P(x) = [I p 0}u(x,0), x(x) = [0 I P }u(x,0). (2.3) 

From (fTTTJ), ((221), and Q one easily gets 

(3(0) = [Ip 0], /3/3* = / P , /3'/3* = = (2.4) 

/3 X * = 0, x(0) = [0 J p ], XV^O, (2.5) 

^(x)=/?'(x) X (x)*. (2.6) 

Put now 

W(A) = {W^(A)}J J=1 = n(/,A)*, (2.7) 

where Wjj are p x p blocks of W. According to (|2.1|) . and (|2.7|1 we 

have [T3]: 

W(A)jW(A)* < j, W(A)*jW(A) < j for 3A < —M. (2.8) 

A pair of p x p meromorphic in the semi-plane 3A < — M matrix functions 
R and Q is called non-singular with j-property if 



R(X)*R(X) + Q(X)*Q(X)>0, [R(X)* Q(xy]j 



R(X) 
Q(A) 



< 0. (2.9) 
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Definition 2.1 liffi Let system M.l)) be given on the interval [0, I] and sat- 
isfy \2.1)) . Then the linear-fractional transformations 

<p{X) = (W 11 (X)R(X)+W 12 {X)Q(X))(W 21 (X)R(X)+W 22 {X)Q(X))~ 1 , (2.10) 

where SA < — M and R, Q are non-singular pairs with j -property, are called 
Weyl functions of this system. 

We denote the set of Weyl functions by J\f(t). From (I2~£1)- (I2~TUJ) it follows 
that 

^(A)V(A) < I p . (2.11) 

The procedure to recover a system from its Weyl functions can be given 
in terms of the S'-nodes introduced in j20]-j2Sl (see also references therein). 
Namely, operator 5* = S(l) acting in the space L 2 (0, /) of squarely integrable 
vector functions and block operator II = [$i $ 2 ], where {k = 1,2) act 
from C p into £-(0, I), are defined in and satisfy operator identity 

AS-SA* = iUU*, A e {Lp(0, 1), Lp(0, 1)}, A = i ■ dt. (2.12) 

Jo 

The construction of these operators is based on the relations: 

s = vz 1 (y^y, $ fc = vz 1 p k , k : = i/3(x) J p(ty ■ dt = V-Av:\ 

(2.13) 

where V- is a triangular and bounded in Lp(0, I) together with its inverse 
operator, are p x p blocks of p. (Here "V- is triangular" means that 
P X V~ = P X V_P X , < x < I, where P x is the orthogonal projector from 
Lp(0, I) onto l/p(0, x).) Moreover operator V~ takes functions with bounded 

derivatives into functions with bounded derivatives, {vi 1 /j(0) = /(0), and 
V- is normalized so that <3>i proves a natural embedding: $ig = g = const. 
Operator $ 2 given in the second relation in ()2.13|) can be presented as a 
multiplication by matrix function and we denote this matrix function by s: 
&29 = s(x)g. It proves [IS] that the following asymptotics takes place for 
the Weyl functions of system (jl.l|) on [0, /]: 

<p(\) = -2i\[ e- 2tXx s(x)dx + 0(\\e~ 2iXl \). (2.14) 
Jo 
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Using (J2.14|) apxp matrix function s(x) with the entries from L 2 (0,l) is 
recovered in ^3] via the Fourier transform 

• poo 

e-"* 8 (x) = — e^A-V(A/2K (A = £ - iv, V > 2M), (2.15) 
where s does not depend on the choice of rj > 2M. 

Theorem 2.2 liffi Suppose if is a Weyl function of system M.l)) on [0, /]. 

T/ien system M.l)) is uniquely recovered from ip by the following procedure. 
First use i2.15\) to introduce operators G {C p , Lp(0,Z)} {k — 1,2) : 

^iS 1 = 9 = const, = s(x)c?, (2.16) 

and operator S = S(l) acting in 1/^(0, 1): 

(2.17) 

These operators are well-defined and bounded. Moreover, S is positive, i.e., 
S > 0, and S is boundedly invertible. 
Next we recover (3* (3 by the formula 

P(x)*P(x) = (irS^r^Il) , (2.18) 

where II is a block operator II := [$i $2], Px is the orthogonal projector 
from Lp(0,l) onto Lp(0,x), and operator S(x) G {Lp(0,x), Lp(0,x)} is given 
by the equality S(x) = P x S(l)P x . 

Finally, the potential v is uniquely recovered from (3*f3 using relations 

Notice that P X U in ()2.18|) is considered as a matrix function [I p s(t)} and 
S(x)~ 1 is applied to this matrix function columnwise. 

There is a simpler way to recover \ from tp so that (3 and v could be 
recovered after that from \. 

Proposition 2.3 Suppose ip is a Weyl function of system M.l\) with a bounded 
on [0, /] potential v. Then matrix function \ defined by the second relation 
in \2.3\) is given by the equality 

X (x) = [0 / p ]-^(,S(x)-V(t))V P s(t)]dt. (2.19) 
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Proof. Denote the upper bound of \ \v\ \ by M, i.e., assume that (|2.1|) holds. 
As ^y_T 1 /J(0) = /(0) according to the first relation in (|2.3|) and second 
relation in (|2.13J) we have: 



s(+o) = (yr^ 2 )(o) = /3 2 (o) = o. 



(2.20) 



Recall that /3 has a bounded derivative and so s = VI 1 fa has a bounded 
derivative too. Apply now operators on the both sides of the third relation 
in 1)2.13)1 columnwise to the matrix function Vs'(x). One easily gets 



Vs(x) = f3(x) I (3{ty(y_s^{t)dt. 



From the second relation in (j2.13j) and from ()2.21|1 it follows that 



(3{x) 



( 


" " 


-f 


Ip 




Jv. 


Jo 


. sity _ 



{yiP x Vs') (t)dt ) = o. 



(2.21) 



(2.22) 



As according to the first relation in ()2.13)) we have S"(x) _1 = V*P X V- formula 
()2.22|) implies /3(x)x(x)* = 0, where £ denotes expression on the right-hand 
side of f!2.19|) . Using S"(x) _1 = V*P X V_ we can rewrite the right-hand side of 
(|2~TTHl as 



[0 I p 



vs')(ty(3(t)dt. 



Thus we have x'(x) = — ^V_s'j (x)*j3(x), which taking into account fix* 



implies x'x* 
equalities 

(3(x)x(xy = 
with the equalities ()2.5|) to see that x 



0. It is immediate also that x(0) = [0 I p ]. Compare 



2(0) = [0 Ip], x\x)x{x) 



L p\-> 
% i-e. 



Tm is true. 



For the case p = 1 formula (J2.19J) have been announced in J^l- I n that case 
(3 is recovered from x — [Xi X2] i n the easiest way: (5 = [x 2 — Xi]- 

Notice also that as matrix function s(x) has a bounded derivative and 
s(0) = 0, formula ()2.17|1 can be written down in a simpler way: 



Sf = f(x) + 



l rx+t 



x—t\ 




r + x — t 



r + t — x 



drf(t)dt. (2.23) 



Finally, let us formulate Borg-Marchenko type theorem. 
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Theorem 2.4 Let ipi and <p 2 be the Weyl functions of the two Dime type 
systems M.l\) on [0, /] with the bounded potentials denoted by V\ and v 2 , 
respectively. Suppose that on some ray SA = c9ftA < and for some < r < 
21 we have 

\\<p 1 (\)- V2 (m=0(\e- ar ) (|A|-oo). (2.24) 
Then V\(x) = v 2 (x) for < x < |. 

Proof. From d2HU) it follows that for each r < 21, |A| — > oo, andSA/|3?A| < 
— 5 < we have 

| MA) - ^(A) 1 1 = -<A jf e~ M ( Sl (|) - s 2 (|)) dt + o(Ae-^)> (2-25) 

where s = si and s = s 2 correspond to systems with v = v% and v = v 2 , 
respectively. Formulas (|2.24j) and ()2.25|) imply the equality 

F(A) ■■= [ z~ lX{t ~ T) -s 2 { t -))dt = 0(l) (2.26) 

on the ray SA = c3?A. Notice that F(X) is bounded in the closed upper 
semi-plane 3A > 0. Using now Phragmen-Lindelof theorem for an angle, in 
view of (J2.26|) we derive that F(X) is bounded also in C_ and thus in the 
whole plane. Moreover, F(X) — ► on the rays in C + , i.e., F(X) = 0. It is 
immediate that Si(x) = s 2 (x) (x < r/2). Applying the procedure to solve 
the inverse problem as in Theorem 12.21 (one can use also formulas (|2.19|) , 
(|2.23Jl ) we finally get v 1 (x) = v 2 (x) for x < r/2. ■ 

Consider now system (jl.ljl on the semi- axis [0, oo). 

Definition 2.5 Let system M.l)) be given on the semi-axis [0, oo). Then 
p x p matrix function <p(X) analytic in some semi-plane 3A < — M is called 
Weyl function of this system if inequalities 

POO 

/ [ <p(\)* I P ]U(X,\)*U(X,\) 

Jo 

hold for all A in the semi-plane. 



dx < oo 



(2.27) 
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Proposition 2.6 \15f System il.lp with bounded on the semi-axis [0, oo) 
potential v 

\\v(x)\\<M (0<s<oo) (2.28) 

has a unique Weyl function (p. Moreover, the matrix disks A/"(/) of Weyl 
functions on the intervals converge to a point, i.e., to an unique function and 
this function proves to be the Weyl function on the semi-axis (p = 0i <oo -N'(l). 

From Theorem 12.41 and Proposition 12.61 follows 

Corollary 2.7 Let ipi and tp 2 be the Weyl functions of the two Dirac type 
systems M-l)) on [0, oo) with the bounded potentials denoted by v\ and v 2 , 
respectively. Suppose that on some ray 3A = c9?A < and for some < r 
equality \2. 2J$ is true. Then t>i(x) = v 2 (x) for < x < |. 

3 Discrete case 

In this section we shall consider skew-self-adjoint matrix discrete Dirac type 
system on the interval: 

W k+1 (X)-W k (X) = ~C k W k (X), C k = C* k = C^\ 0<k<n. (3.1) 

When p = 1, then either C k = ±I 2 or C k = U(k)jU(k)* (U(k)U(k)* = 
I2). If C k = U(k)jU(k)* (k > 0), system (|3.1|) is auxiliary system for 
isotropic Heisenberg magnet model. Therefore here we also assume that 
C k = U(k)jU(k)* and matrices U(k) are unitary, i.e., 

C k = I 2p -2(3(k)*(3(k), (3(k)(3(ky = I p , 0<k<n, (3.2) 

where (3{k) = [ftx{k) P 2 (k)] = [0 I p ]U(k)* are p x 2p matrices with p x p 
blocks /3i(k), f3 2 (k). Introduce now simple additional conditions: 

det /3i(0) ^ 0, det (3{k - l)/3(k)* ^0, < k < n. (3.3) 

Similar to the continuous case we shall define Weyl functions of the system 
via Mobius (linear-fractional) transformation 

<p{\) = (W n (A)i?(A) + W 12 (A)Q(A)) (W 21 (A)i?(A) + W 22 (A)g(A))" 1 , (3.4) 
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where we put 

W(A) = {VVy(A)}?j =1 = W n+1 (A)*. (3.5) 

Here 2p x 2p solution of (|3.1jl is normalized by the condition Wo (A) = / 2 p 
and coefficients of the Mobius transformation are p x p blocks of W. 
We shall be interested in the properties of (p(X) in the neighborhood of A = 
i. So we require that R and Q are p x p matrix functions analytic in the 
neighborhood of A = % and such that 

det (w 21 {i)R{i) + W 22 (i)Q(z)) ^ 0. (3.6) 

Such pairs R, Q always exist as the rows of [W 2 i(z) VV 2 2(z)] are linearly 
independent: 

rank[W 21 (*) W 22 (i)] = p. (3.7) 

Indeed, to prove ()3.7|) we can take into account (j3.1j) and Ck = U(k)jU(k)* 
and derive 

n 

W(-i) = 2 n+1 H x (k)*x(k), x (k) = [ Xl (k) X 2(k)] = [I p 0]U{k)\ 

k=0 

(3.8) 

Notice that by U(k)*U(k) = I 2p we have 

X (k)P(k)* = 0, X (k)x(k)* = I P . (3.9) 
From (jHS) and (HTEJ) it follows that 

[W 2 i(z) W 22 (z)] = 



= 2™ +1 X2 (0)*(x(0)x(l)*) (x(l)x(2)*) • • • (x(« - l)x(n)*)x(n). (3.10) 

It remains to show that inequalities ()3.3|) imply: 

detx2(0)^0, det - l)x(&) V °> < A; < n. (3.11) 

Suppose det X2(0) = 0. Then there is a vector / ^ such that /*x 2 (0) = 
0. So, in view of the first relation in (j3.9j) we get = /*x(0)/3(0)* = 



/* Xi(0)A(0)*. As det /3i(0) ^ it follows that /* Xi(0) = and thus 
/*x(0) — 0. But according to the second relation in (|3.9|) the lines of x(0) 
are linearly independent and we come to a contradiction. The first inequality 
in (EHTJ) follows. 

Suppose there is a vector / ^ such that f*x(k ~~ — 0- By ()3.9|) 

we then obtain 

r x (k-i) = rm, /Vo. 

Hence, taking into account the second inequality in (J3.3j) we have 

f* X (k - l)(3(k - 1)* = f*/3(k)/3(k - 1)* ^ 0. 

As the first equality in (J3.9)) yields x{k — l)ft(k— 1)* = 0, we get contradiction, 
i.e., the second inequality in (j3.1H) is valid too. Now relations (j3.10J) and 
(l3~TTJ) imply equality (l3~7|) . 

Definition 3.1 Let system \3. 1)) be given on the interval < k < n and 
satisfy \3.$fy . \3. 3]) . Suppose R and Q are p x p matrix functions analytic in 
the neighborhood of A = i and such that inequality A3.6)) holds. Then linear- 
fractional transformations ip of the form \3.4\ l are called Weyl functions of 
this system. The pair R, Q satisfying our conditions is called admissible. 

Example 3.2 Put n — 1, Cq — —j, C\ = J, where 



J 



I 



KjK* 



K 



1 

7! 



K* = K 



-i 



(3.12) 



In view of A3.1ty) we have Cq = JjJ, and so U(0) = J. By the second equality 
in KS.lty we get U(l) = K. R follows that 



/?(o) = [o i p ]u(oy = [i p o], = [o i p p(iy = -y=[-i p i p \. 



Thus the conditions A3, ty) and \3. 3\) are fulfilled. Moreover we have 



(3.13) 



W 2 {\) 



•2p 



i.e., 



W(X)=l2 P -j(j-J)-^Jj. 



(3.14) 
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In particular, we have W^i^) = 2, = 2. Thus condition \3. 6)) tafces 

t/ie /orm 

det (fl(i) + Q(i)) ^ 0. (3.15) 

( | J. i^| ) Mobius transformation \3.J$ , that describes Weyl functions, takes 
the form 

ip(\) = ^ zl (\R + iQ)(iR + \QT l . (3-16) 
A + % 

Denote the set of Weyl functions of system (|3.1|) on the interval < k < n 
by M(n). 

Remark 3.3 The sets of Weyl functions introduced by Definition IJ.il are 
decreasing. That is, given system \3. 1\) we have M(l) ~3 M{n) for I < n. 
Indeed, the pair 



R(X) 
Q(A) 

satisfies condition 



hp + J C i+t ) x • • • x ( hp + ~C, 



R{\) 



(3.17) 



w l+1 (\y 



R(X) 
Q(A) 



W(A) 



R(X) 

L QW J 



(3.18) 



So if the pair R, Q is admissible for system \3. 1\) on the interval < k < n, 
then the pair R, Q is admissible for system \3.1\) on the interval < k < I. 
Moreover Weyl function tp of system \3.1\) on the interval < k < n that is 
determined by R and Q coincides with the Weyl function of system \3. 1)) on 
the interval < k < I that is determined by R and Q. 

The next theorem solves inverse problem to recover system (|3.1|) from its 
Weyl function. 

Theorem 3.4 Suppose <p is a Weyl function of system \3.1\) satisfying con- 
ditions \3. and \3. 3\) . Then system \3. 1)) is uniquely recovered from the first 

■1 + z* 



n+1 Taylor coefficients {ak}% =0 of ip 
procedure. 



at z = by the following 
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First introduce (n + l)p x p matrices $i ; $2 : 









$1 = 


h 


, $2 = - 









Q'o 

a + «i + • • • + ct n 



(3.19) 



JTien introduce (n + l)p x 2p matrix IT and (n + l)p x (n + l)p lower 
triangular matrix A by their blocks : II = [$1 $ 2 ], 



A:=A(n) = {aj-*}*j= ' ttr 




(3.20) 



iVext we recover (n + l)p x (n + l)p matrix S as a unique solution of the 
matrix identity 

AS - SA* = iim*. (3.21) 

This solution is invertible and positive, i.e., S > 0. Finally matrices /3(h)* /3(h) 
are easily recovered from the formula 



ITS'- 1 !! = B*B, B := Bin) 



/3(0) 

m 

(3(n) 



(3.22) 



Now matrices Ck and system hS. 1\) are defined via hS.tyi . 

Proof. Step 1. The method of the proof coincides with the method of the 
proof of Theorem 12.21 Put 



K{r) 



K (r) 
K x (r) 

KJr) 



(3.23) 



where Kj(r) are p x (r + l)p matrices of the form 

K j (r) = ip(j)WY — M-iY PUT/ 2 o.-.o]- 



(3.24) 
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From (jnUgP - flPS) it follows that 

K(r) -K{r)* = iB(r)B(r)*. (3.25) 

By induction we shall show in the next step that K is similar to A: 

K(r) = y_(r)A(r)y_(r) _1 (0 < r < n), (3.26) 

where VL(r) are block lower triangular matrices. Taking into account 

()3.26|) and multiplying both sides of ()3.25|) by y_(r) _1 from the left and 
by (V_(r)*) from the right we get 

A{r)S{r) - S{r)A{r)* = iU(r)U(r)*, (3.27) 



S(r) :=y_(r) _1 (VL(r)*) \ n(r) := y_(r) -1 S(r). (3.28) 

Moreover, Step 3 will show that matrix VL(n) can be chosen so that the 
equality 

n = [$! $ 2 ] = V-in^Bip) (3.29) 

holds, i.e., II = U(n). (Here $i and $2 are given by (|3.19|l .) 

Identities f!3.27|) have unique solutions S{r) as the spectra of A and A* 
don't intersect. In particular, by (j3.21|) and (|3.27|) one can see that S := S{n). 
Hence we derive from ()3.28j) and ()3.29j) that S > and the first equality in 
(I3~22J) holds. It remains only to prove (J3~2T)|) and (J3~2U1) . 

Step 2. Now we shall consider block lower triangular matrices V-(k) 
(0 < k < n): 



K_(0) = «_(0) = A(0), \/_(fc) 



X(fc) 







(k > 0), (3.30) 



where f_(/) are p x p matrices, X(k) = [X (k) X(k)] is p x kp matrix, 
Xo(k) is an arbitrary pxp block, and X(k), V-(k) are given by the formulas 



X{k) =i [p(k)[P(0y ...p(k- l)*]V-(k 



'(fc-i)p 




- v-(k)[I p ...I p ] 



x(A(A;-2) + ^/ (fc _ 1)p ) \ v-(k) = (3(k)(3(k - l)*v-(k - 1). (3.31) 
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According to (|3.2(Jj) we have A(0) = (i/2)I p . From the second relation in 
(jSHD and definitions (TT23|) and (jHUlD it is immediate that K(0) = {i/2)I p 
and so ()3.26|) is valid for r = 0. Assume that ()3.26|) is true for r = k — 1, and 
let us show that ()3.26|) is true for r = k too. It's easy to see that 



V-(k) 



-1 



v^{k-iy l o 

-v_(A;) -1 X(A;)y_(A; — v^k)' 1 



(3.32) 



Then in view of definitions ()3.20|) and ()3.30|) our assumption implies 



V-(k)A(k)V-.(k)~ 1 



K(k-l) 
Y(k) 



-I 



(3.33) 



where Y(k)= (X(k)A(k - 1) + iv-(k)[I p . . . I p }) ^v-(k) 

V^k-l)- 1 
X _ —v-(k)~ 1 X(k)V-.(k - l) -1 _ ' 

Rewrite product on the right hand side of the last formula as 

Y(k) = (x(k)(A(k - 1) - % -I kp ) + iv-(k)[I p ...I p ])V4k - (3.34) 
From (13~2TH) and (EPIJ) it follows that 



Y(k) = ^X{k)(A{k-2) + -I (k „ l)p )+iv„{k)[I p . . .L p }} iv_{k)\V-(k-iy\ 

(3.35) 

Notice that the sequence [I p . . . I p ] of identity matrices in (|3.35|) is one block 
smaller than in ()3.34|) . By (|3.31|) and ()3.35|1 we have 



Y(k) = i/3(k) \p{p)*...p(k-l)*\V-(k-l) 



! (fc-l)p 





P(k-l)*v-(k-l) 



xV^k-iy 1 . 
Finally formulas ()3.3()|1 and ()3.36j) imply 

Y(k) = i(3(k)[P(0)* . . . (3(k - 1) 



(3.36) 



(3.37) 
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According to the second relation in (J3.2)) and formulas ()3.24)) and (J3.37)) we 

get 



Y(k) -I P 



K k {k). 



(3.38) 



Using now (J3.23)) and (|3.37j) one can see that the right hand side of (|3.33j) 
equals K(k). Thus 1)3.26)1 is true for r = k and therefore for all < r < n. 

Step 3. To derive ()3.29)) we shall first prove that matrices VL(r) given by 
(13.30)1 and (j3.31j) can be chosen so that 



V^rY'B^r) 



Bx{r) := B{r) 







Pi(0) 
L Pi(r) J 



(3.39) 



In other words the arbitrary till now blocks X (r) can be chosen so. Indeed, 
by definition in (j3.22j) and the first equality in (J3.30)) formula ()3.39)) is true 
for r = 0. Assume that ()3.39)) is true for r = k — 1. Then from ()3.32j) it 
follows that ()3.39j) is true for r = k if 



-v^(k)~ 1 X(k) 



L Ip J 



+ v.(k)- 1 p 1 (k) = I p . 



(3.40) 



It implies that we get equality ()3.39|) for r = k putting 



X (k)=/3 1 (k)-v4k)-X{k) 



(3.41) 



Hence by a proper choice of matrices Xo(r) we obtain ()3.39)) for all r < n. 
It remains to prove that 



V4n)- 1 B 2 {n) = $ 2 , B 2 (n) 



L ft(n) 



(3.42) 



For that purpose we shall consider matrix function W n+1 (A), that is used 
in (13.5)) to define coefficients of the Mobius transformation (j3.4j) (and thus 
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Weyl functions). Namely we shall prove so called transfer matrix function 
representation of W n+ i(X): 

where 

w A (r, A) = I 2p - tUirySiry^Air) - A/ (r+1)p ) _1 n(r). (3.44) 

Transfer matrix functions of the form 1)3.44)1 have been introduced and studied 
by Lev Sakhnovich [2D]-[23]- In particular, identity ()3.27J1 implies 

w A (r,n)*w A (r, A) = I 2p 

X)U(ry(A(r)* -flI (r+1)p )- 1 S(r)- 1 (A(r) - Xl^y^r). (3.45) 

Moreover, according to factorization Theorem 4 from [20] (see also [221, P- 
188) we have 

w A (r, A) = (l 2p -iIL(r)*S(r)- 1 P*(PA(r)P*- XI p )~\PS(r)- 1 P*y 1 

xP5(r) -1 n(r)j?/;4r-l,A), P = [0 ... I p }. (3.46) 
Taking into account (|3.2U|) . ()3.28|) . and ()3.30|) we obtain 

(PA(r)P* - M.y 1 = C-- A)- 1 ^, P5(r)- x P* = u_(r)*^(r), (3.47) 

PSir)' 1 ^) = v_(r)*PB(r) = u_(r)*/3(r). (3.48) 
Substitute (f3~¥7f) and (EP%|) into lEOE|) to get 

«A(r, ~) = (/ 2p - -^/5(r)*/?(r)) WA (r - 1, ^). (3.49) 
From the definitions (|3.2U|) . (|3.28|) . and ()3.44|) we also easily derive 

w A (0, = J 2p - JLB(0)*B(0) = h P - t^-/3(0)*/3(0). (3.50) 
2 z — A « — A 
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On the other hand system (|3.1|) with additional condition (|3.2j) can be rewrit- 
ten as 



W{r + l,X) 



X-i 



2i 



A V i — A 

In view of the normalization W(0) = / 2p formulas ()3.49|) - ()3.51|) imply ()3.43|) . 
From (g3H)) and (|3~4"5l) it follows that 



P(r)*/3(r))W(r,\). (3.51) 



W(n + l,X)W(n+l,\y 



■ \ n+l / \ i ■ \ n+1 

A — t\ I X + 1 



X J \ A 

Let us include now Weyl functions into consideration and put 



A{\) := 



A 



X-i 



2n+2 



[y?(A)* J p ]W(n + l,A)"W(n + l,A) 



According to (|3.4|) . ()3.5|) . and ()3.52|) we have 



-4(A) 



A + i 



A 



2n+2 



((W 21 (A)i?(A) + W 22 (A)Q(A))*) 



*\ -i 



x^(A)* J R(A) + g(A)*g(A)J(W 21 (A) J R(A) + W 22 (A)g(A)) \ (3.54) 
By ()3.6|) and ()3.54|) A is bounded in the neighborhood of A = i: 

\\A{X)\\=0{1) for A (3.55) 

Substitute now ()3.43|) and ()3.45|) into (j3.53|) to obtain 

.4(A) = [<£>(A)* I,](l 2 , + A)n(n)*(i(«)' - ^hn^Y'Sin)- 1 



(3.52) 



(3.53) 



A 



x { A ( n ) - ^ J (n+i) P ) ln (^) 



<p(\) 



(3.56) 



Notice that 5"(n) > 0. Hence formulas ()3.55|) and ()3.56|) imply that 



¥>(A) 



0(1) for A -> i 



(3.57) 
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Recall that U(n) = V-(n) l B(n) and A(n) is denoted by A. Represent now 
II(n) in the block form 

II(n) = [$i(n) $ 2 (n)], $*(n) = VLfa)- 1 ^) (A; = 1,2). (3.58) 

According to (|3.19j) and (J3.39j) we have $i(n) = $x- Hence multiplying the 

matrix function on the left hand side of ()3.57j) by (A—~I( n+1 ) p ) j $* 
we derive 



A^z. 



(3.59) 



The matrix A— \l< n +i)p is easily inverted explicitly (see, for instance, formula 
(1.10) in J7j). As a result one obtains 



A 



l (n+l)p 



-1 



n „n— 1 



z-A 



[q n 4 



q U <? : =^f ( 3 - 60 ) 



Moreover we get 



2 V +1 -/ P )fg-/ P ^' 



(3.61) 



Putting A = i[ — \ i.e., z = ( \ '■ ) we derive from ()3.61|) that 



V-zJ VA • / 

A • 1 \-i 



(^(A-^Dp)- 1 ^)" =(-z^ +1 + 0(z 2 " +2 ))J p (z-0). (3.62) 

Taking into account ()3.60|) and ()3.62|) we rewrite ()3.59|) as 
1 + z- 



<p\i 



l-z 



(l-z)[I p zl p z 2 I p ...]$ 2 (n) 



0(z n+1 ) (3.63) 



for * -> 0. From (j3~TOj) and (f3~K3~j) follows that $ 2 (n) = $ 2 , i.e., flgZj) is 
true. As = $i and $ 2 (n) = $ 2 , so n(n) = II and formula ()3.29|1 is 

finally proved. ■ 



From Theorem 13.41 and Remark 13.31 follows Borg-Marchenko type result. 
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Theorem 3.5 Let <p and (p be Weyl functions of the two discrete Dirac type 
systems \3.1\) that satisfy conditions \3. ty) and \3. 3\) . Denote by Ck (0 < 
k < n) the potentials Ck of the first system and by Ck (0 < k < n) the 
potentials of the second system. Denote Taylor coefficients of cplilj^-)) 



and <f\i(jz^JJ at z = by and {o^}, respectively, and assume that 



a.k = 0(k for k < I (I <n, I <n). Then we have Ck = Ck for k < I. 

Proof. According to Remark 13.31 tp and (p are Weyl functions of the first 
and second systems, respectively, on the interval < k < I. By Theorem 13.41 
these systems on the interval < k < I are uniquely recovered by the first 
I + 1 Taylor coefficients of the Weyl functions. ■ 

Step 3 of the proof of Theorem 13.41 implies the following corollary. 

Corollary 3.6 Weyl functions of system \3.1\) that satisfies conditions \3.ty) 
and \3. 3\) admit Taylor representation 



where ipk ar & P x P blocks of $2(^) = {^aJaUo- Here $2(^) is given by the 
second equality in \3.5§) . where V-(n) is defined by formulas \3. 3U\) . \3. 31]) . 
and (3411 ). 

Moreover, from the proof of Theorem 13.41 follows a complete description of 
the Weyl functions in terms of the Taylor coefficients. 

Theorem 3.7 (i) Let system \3. 1)) be given on the interval < k < n 
and satisfy \3. fy) . \3. 3\) . Then analytic at A = i matrix function ip is a 
Weyl function of this system if and only if it admits expansion \3. 64\ ), where 
matrices ipk are defined in Corollaru \3. 61 

(ii) Suppose ip is a p x p matrix function analytic at X = i. Then ip is a 
Weyl function of some system \3. 1)) given on the interval < k < n and 
satisfying \3. \3. 3\) if and only if the matrix S uniquely defined by the 
identity AS — SA* = iUU* is invertible. Here IT = [$i $2] is given by 
\3.19\) . where {a^} are Taylor coefficients of ipiii ) ) at z = 0. 






(i (j^-) )=-^ + ^0-^l)z + ... + (Vn-l - ^n)z n + 0(z n+1 ) (* - 0), 

(3.64) 
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Proof. Suppose tp is apxp matrix function analytical at A = i and satisfying 
(I3.64j) . and put 



R(X) 



:= W(A)- 1 



tp(X) 



(3.65) 



Fix a Weyl function ip of system ()3.1|) . and denote by i?, Q some admissible 
pair that grants representation (J3.4j) of Rewrite (|3.65|) in the form 



R(\) 
Q(A) 



W(A) 



-i 



£(A) 



+ W(A) 



-i 



¥>(A) - £(A) 




(3.66) 



As ip is the Mobius transformation ()3.4|) of the admissible pair R, Q, we have 



W(A) 



-i 













_ 0(A) _ 



(W 21 (A)i?(A) + W 22 (A)Q(A))" 



(3.67) 



Thus the first summand on the right-hand side of (|3.66j) is analytic at A = 
i. Taking into account that expansion ()3.64|) is valid for ip and tp we de- 
rive ip(\) - (p(X) = 0((A - i) n+1 ) for A i. From (g32I) it follows also 

that W(A) _1 = — r — — — r — — Win + 1,A). So the second summand 

v ' (A + i) n+1 ( \ - i ) n+1 y ' 

on the right-hand side of ()3.66|) is analytic at A = i too. Therefore the 
pair R, Q is analytic at A = i. Moreover, according to ()3.65j) we have 
W 2 i(A)i?(A) + W 22 (A)Q(A) = I p . Hence, inequality (l3~T)J) holds and the pair 
R, Q is admissible. It easily follows from (|3.65j) that <p admits representation 
()3.4|) with this R, Q, i.e., <p is a Weyl function of our system. 

Vice versa, if ip is a Weyl function of our system, then by Corollary 13.61 
the expansion (|3.64j) is true. The statement (i) is proved. 

According to the proof of Theorem 13.41 when tp is a Weyl function, then 
matrix S uniquely defined by the identity AS — SA* = iUU* is invertible. It 
remains to show, that if S is invertible, then <p is a Weyl function. For this 
purpose introduce notations 



Pr= [0 



I p ], P r = [I (r+1)p 0], S(r) 



P r S P r , 



nw 



p r n, 

(3.68) 
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where P r are p x (r + l)p and P r are (r + l)p x (n + l)p matrices. Assume 
that det S 7^ 0. As the identity (|3.21|) is equivalent to 

5(A*-A/ (n+1)p )- 1 -(A-A/ {n+1)p )- 1 5 = z(A-A/ (n+1)p )- 1 nn*(A*-A/ (n+1)p )- 1 , 

using residues we standardly derive 
1 f°° 

S = tt ( A ~ A/(„ + i )p )- 1 nn*(A* - Xl^y^dX. (3.69) 

^ J -co 

Therefore from det S ^ it follows that S > 0, and so we have S*(r) > 
and det S(r) ^ 0. Put 

p(r) := (P r 5(r) _1 P r ) _ 5p ri 5( r )-i n ( r ) (0 < r < n). (3.70) 

Matrices /3(r) satisfy conditions ()3.2|) and ()3.3|) . Indeed, from ()3.7()j) we get 

/3{r)/3{r)* = (P r 5(r)- 1 P r )"^P r 5(r)- 1 n(r)n(r)*5(r)- 1 P;(P r 5(r)" 1 P r )"i 

(3.71) 

From (EHU) follows identity (jS23), i.e., 

S , (r)- 1 n(r)n(r)*S'(r)- 1 = -z(5(r)^M(r) - A(r)*5(r)- X ). 

It is also true that P r A(r)* = — ~P r and A(r)P* = |P r . Hence we can rewrite 
(I3~TTJ) as 

p{r)p{r)* = -{(PrSir^Pry^Prfsir^Air) - A(r)*S(r)- x ) 

xP^PrSiry'Pr)^ =I P . (3.72) 

Thus the condition from ()3.2|) on /? is proved. Notice further that according 
to (l3~m and (l3~70l we have 

A(0) = (S-CO)- 1 )"^^)- 1 = S(0)-$, (3.73) 

and so the first inequality in ()3.3|) is also true. To prove the second inequality 
in ()3.3j) we should write down 5(r) _1 (r > 0) in the block form: 

( _ [ £(r - I)" 1 + 5(r - l)- 1 S 12 tS 21 S(r - l)" 1 -5(r - iy 1 S 12 t ' 
[T) [ -tS 21 S(r - l)- 1 tj' 

(3.74) 
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where t = (s — S 2 iS(r — 1) 1 S\ 2 ) , and Si 2 , S 2 \ and s are blocks of S(r): 



S(r) 



S(r-1) S 12 
S21 s 



(3.75) 



Formula ()3.74|) is easily checked directly. (Conditions det S(r) 7^ and 
det S(r — 1) 7^ imply invertibility of s — S 2 xS(r — l) _1 S'i 2 ; see, for instance, 
j2Sj, P- 21.) In view of (|3.7U|) and (|3.74|) the second inequality in (|3.3|) is 
equivalent to 

det ([-S 21 S{r - I)- 1 I p ]U{r)U{r - l)*S{r - l^P*^) ^0 (r > 0). 

(3.76) 

Recall that S(r — 1) also satisfies operator identity: 

A(r - l)S(r - 1) - S(r - l)A(r - 1)* = iU(r - l)U(r - 1)*. 
Therefore similar to the proof of ()3.72|) one gets 



[S 21 S(r - l)- 1 I p ]U(r)U(r - l)*S(r - l^P^ = ~S 21 S(r - l)"^ 



-iS 21 A(r - l)*S(r - l)" 1 ^ + [I p W]n(r ~ l)*3(r - l)"^, (3.77) 

where {ipk} are the blocks of $2- Consider now the first r blocks in the lower 
block rows on the both sides of the identity (|3.27|) . In view of (J3.75|) we have 

l -S 21 +i[I p ... I p ]S(r- 1) - S 21 A(r - 1)* = i[I p ^ r ]U(r - 1)*, i.e., 



i(^S 21 -S 21 A(r-l)*-i[I p VrMr-l)*)^-!)-^;.^/,. (3.78) 

By (jTZID and (jTZj) inequality (l3~7oT) holds. 

So formula ()3.70|) defines system (j3.1|) that satisfies conditions (|3.2J) and 
()3.3j) . Moreover, by Theorem 4 [20] the transfer matrix functions ^(r, A) 



22 



corresponding to our matrices S(r) admit factorizations (J3.46|) . In view of 
(I3.7U|) this implies 



\ 9 ' 

Mn, ? ) = II (> - T^rftrypir) 



r=0 



i - A 



(3.79) 



Compare (|3.79|) with (|3.51|) to derive for the fundamental solution of the 

/ ,\n+l 

constructed system the equality W(n + 1, A) = ( ^ j ^(n, |), i.e., 



/ \ \ 71+1 \ 

W(A)-=( — ) «*(»,-). 



A, 
2' 



(3.80) 



Let us consider now y(A) and prove ()3.57|) . It is easy to see (|17j. p. 452) 
that the fc-th block row of (A — f i( n +i)p) is given by the formula 

i(l-z) 



T k (X) 



yk + l 



X 



[l-z)I p {l-z)zl p 



[1 - z)z k - 1 I p z k I p 



where z 



X-i 
X + i 



. It follows also that 



(A - -/ (n +l)p) $1 = ~ COI 



± p a, ± p 



z~ n L 



where col means column. By ()3.8H) and ()3.82|) we have 



ip(\) 



(3.81) 



(3.82) 



%{\-z)( fjl + z 



fT^)) + + (i>i - ^o)z + ... + (if> k - ^ k -i)z h 



k=0 



(3.83) 



Formulas (J3~64jl and (J3~83jl imply flEgfl . From (jTHjl . (pTBTj) . and (J3~80jl it 

¥>(A) 



follows that 



W(A)- 1 



0(1) for A — > «. Thus the pair R, Q given 



by ()3.65|) is admissible. As (p admits representation ()3.4|) with R, Q given by 
()3.65|) . so ip is a Weyl function of the constructed system. ■ 
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Example 3.8 According to A3. 16]) for system considered in Example 



we 



have «o = and a\ = I p . Hence, by Theorem TfT) the set of Weyl functions 
ip of this system is defined by the expansion 



<p\i 



l + z 



z + 0(z 2 



0. 



(3.84) 
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